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Abstract We study the properties of charmed pseudoscalar and vector mesons in 
dense matter within a unitary meson-baryon coupled-channel model which incorporates 
heavy-quark spin symmetry. This is accomplished by extending the SU(3) Weinberg- 
Tomozawa Lagrangian to SU(8) spin- flavor symmetry and implementing a suitable 
flavor symmetry breaking. Several resonances with negative parity are generated dy- 
namically by the s-wave interaction between pseudoscalar and vector meson multiplets 
with l/2 + and 3/2 + baryons. Those states are then compared to experimental data as 
well as theoretical models. Next, Pauli-blocking effects and meson self-energies are in- 
troduced in a self-consistent manner to obtain the open-charm meson spectral functions 
in a dense nuclear environment. We finally discuss the formation of D-mesic nuclei. 
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1 Introduction 

The interest on the properties of open and hidden charmed mesons was triggered in 
the context of relativistic nucleus-nucleus collisions in connection to the charmonium 
suppression [I] as a probe for the formation of Quark-Gluon Plasma (QGP). The 
primary theoretical effort is to understand the interaction between hadrons with the 
charm degree of freedom. The nature of new charmed and strange hadron resonances 
is an active topic of research, with recent data coming from CLEO, Belle, BaBar and 
other experiments [2]. In the upcoming years, the experimental program of the future 
FAIR facility at GSI [3] will move from the light quark sector to the heavy one and 
will face new challenges where charm plays a dominant role. 

Approaches based on coupled-channel dynamics have proven to be very successful 
in describing the existing experimental data. In particular, unitarized coupled-channel 
methods have been applied in the meson-baryon sector with charm content 4,5,6,7, 
8,9 , partially motivated by the parallelism between the yl(1405) and the /l c (2595). 
Other existing coupled-channel approaches are based on the Jiilich meson- exchange 
model |1UII11| or on the hidden gauge formalism [12| . 

However, these models are not fully consistent with heavy-quark spin symmetry 
(HQSS) 13,14,15 , which is a proper QCD symmetry that appears when the quark 
masses, such as the charm mass, become larger than the typical confinement scale. Aim- 
ing at incorporating HQSS, an SU(8) spin-flavor symmetric model has been recently 
developed [16II17] . similarly to the SU(6) approach in the light sector of Refs. [1811 19] . 
The model can generate dynamically resonances with negative parity in all the isospin, 
spin, strange and charm sectors that one can form from an s-wave interaction between 
pseudoscalar and vector meson multiplets with 1/2 -1- and 3/2 + baryons [20] . 

In this work we will focus on the modifications of the dynamically-generated states 
in the nuclear medium and on the properties of open-charm mesons in dense nuclear 
matter. We will study the open-charm meson spectral functions in this nuclear en- 
vironment within a self-consistent approach in coupled channels, and discuss their 
implications on the formation of D-mesic nuclei. 

2 SU(8) Weinberg- Tomozawa interaction with heavy-quark spin symmetry 

In QCD all types of spin interactions vanish for infinitely massive quarks: the dynamics 
is unchanged under arbitrary transformations in the spin of the heavy quark. This is the 
prediction of HQSS. Thus, HQSS connects vector and pseudoscalar mesons containing 
charmed quarks. Chiral symmetry fixes the lowest order interaction between Goldstone 
bosons and other hadrons in a model independent way; this is the Weinberg- Tomozawa 
(WT) interaction. Then, it is very appealing to have a predictive model for four fla- 
vors including all basic hadrons (pseudoscalar and vector mesons, and l/2 + and 3/2 + 
baryons) which reduces to the WT interaction in the sector where Goldstone bosons 
are involved and which incorporates HQSS in the sector where charm quarks partic- 
ipate. Indeed, this is a model assumption which is justified in view of the reasonable 
semiqualitative outcome of the SU(6) extension in the three-flavor sector [21] and on 
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a formal plausibleness on how the SU(4) WT interaction in the charmed pseudoscalar 
meson-baryon sector comes out in the vector-meson exchange picture. 

The model obeys SU(8) spin-flavor symmetry and also HQSS in the sectors where 
the number of c- and c- quarks are conserved separately. This is a model extension of 
the WT SU(3) chiral Lagrangian 20 . Schematically, 

4 U t 8 ' = JI H m1 ® M l63 a ® [B ] ® B]e»]i, (1) 

which represents the interaction between baryons (in the 120 irrep of SU(8)) and mesons 
(in the 63) through t-channel exchange in the 63. In the s-channel, the meson-baryon 
space reduces into four SU(8) irreps, from which two multiplets 120 and 168 are the 
most attractive. As a consequence, dynamically-generated baryon resonances are most 
likely to occur in those sectors, and therefore we will concentrate on states which belong 
to these two representations. 

The SU(8)-extended WT meson-baryon interaction is given by 



\f,Sj V 2M i V 2M j 



Here, i (j) are the outgoing (incoming) meson-baryon channels. The quantities Mi, Ei 
and fi stand for the baryon mass and energy, in the center of mass frame, and the meson 
decay constant in the i channel, respectively. are the matrix elements coming from 
the SU(8) group structure of the coupling for the various charm, strange, isospin and 
spin (CSIJ) sectors, which display exact SU(8) invariance. However, this symmetry 
is severely broken in nature, so we implement a symmetry-breaking mechanism. The 
symmetry breaking pattern follows the chain SU(8) D SU(6) D SU(3) D SU(2), where 
the last group refers to isospin. The symmetry breaking is introduced by means of a 
deformation of the mass and decay constant parameters. This allows us to assign well- 
defined SU(8), SU(6) and SU(3) labels to the resonances and to find HQSS invariant 
states. 

We solve the on-shell Bethe-Salpeter equation in coupled channels so as to calculate 
the scattering amplitudes, Ty , by using the interaction matrix V as kernel: 

T ^ = i-vwc W y W- ® 

G(s) is a diagonal matrix containing the meson-baryon propagator for each channel. 
D, T, V, and G are matrices in coupled-channel space. The loop function G(s) is 
logarithmically ultraviolet divergent and it is regularized by the means of a subtraction 
point prescription [22] . 

The poles of the scattering amplitudes are the dynamically-generated baryon res- 
onances. Close to the pole, the T-matrix behaves as 

Tij (s)K J i9 * , (4) 



with the mass (m/j) and the width (-T^j) given by y^sTj = m^j — irpj/2, while is the 
coupling to a given meson-baryon channel. 
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3 Charmed and strange dynamically-generated baryon states 

Dynamically generated states in different charm and strange sectors are predicted 
within our model 16,17,20 . Some of them can be identified with known states from 
the PDG |23| . This identification is made by comparing the PDG data on these states 
with the mass, width and, most important, the coupling to the meson-baryon channels 
of our dynamically-generated poles. Next we will concentrate on the (7=1, 5 = 
sector since the states generated in this sector are the most relevant ones for the latter 
analysis of open-charm mesons in dense nuclear matter. 

We reproduce the three lowest-lying states of Ref. [16] in the A c sector, which come 
from the most attractive SU(8) representations. However, those states have slightly dif- 
ferent masses due to the different subtraction point, and the use of slightly different D s 
and Dg meson decay constants. The experimental yl c (2595) resonance can be identified 
with the pole that we obtain around 2618.8 MeV, as similarly done in Ref. [16j. The 
width in our case is, however, bigger because we have not changed the subtraction 
point to fit its position as done in Ref. [TB]. Nonetheless, the dominant three-body 
decay channel A c mr [23] is not included in our calculation. A second broad A c reso- 
nance at 2617 MeV is, moreover, observed with a large coupling to the open channel 
S c n, very close to yl c (2595). This is the same two-pole pattern found in the charmless 
/ = 0, S = — 1 sector for the yl(1405)|24j. A third spin- 1/2 A c resonance is seen around 
2828 MeV and cannot be assigned to any experimentally known resonance. With regard 
to spin-3/2 resonances, we find one located at (2666.6 — i26. 7 MeV). As seen in Ref. If). , 
this resonance is assigned to the experimental yl c (2625). The i-channel vector-exchange 
model of Ref. [7] also observed a similar resonance at 2660 MeV. The novelty of our 
calculations with respect to that one is that we obtain a non-negligible contribution 
from the vector meson-baryon channels to the generation of this resonance. 

For C — 1,S — 0,1 — 1,J= 1/2 (S c sector), three E c resonances are obtained 
with masses 2571.5, 2622.7 and 2643.4 MeV and widths 0.8, 188.0 and 87.0 MeV, 
respectively. Those states nicely agree with the three lowest lying resonances found in 
Ref. [16] and are pure predictions of our model. Ref. [9] predicts the existence of two 
resonances in this sector. However, only one of them can be identified to one of ours but 
with a strong vector meson-baryon component. Moreover, we predict two spin-3/2 S c 
resonances. The first one, a bound state at 2568.4 MeV, is thought to be the charmed 
counterpart of the £'(1670). The second state at 2692.9 — i33.5 MeV has not a direct 
experimental comparison. 



4 Open-charm mesons in nuclear matter 

We start this section by first analyzing the behavior of the dynamically-generated 
states within the C = 1 and S = sector in dense nuclear matter. Then, we study the 
properties of open charm (D and D*) mesons in a dense nuclear environment related 
to the modification of dynamically-generated A c and E c states in this environment. 

The self-energy and, hence, spectral function for D and D* mesons are obtained 
self-consistently in a simultaneous manner, as it follows from HQSS by taking, as bare 
interaction, the SU(8)-extended WT interaction described in Sec. [2] We incorporate 
Pauli blocking effects and open charm meson self-energies in the intermediate propa- 
gators [25] for the in-medium solution. 



5 




The D and D* self-energies are obtained by summing the transition amplitude for 
the different isospins over the Fermi sea of nucleons, ra(p): 



n D (qo,q) = J -^y3 n(p) [T DN 
d 3 p 



(7=0,J=l/2) 



3T DN 



(I=l,J=l/2), 



2 (7=0,J=3/2) 

3 D * N 



1 „(I=0,J=l/2) ^ T (I=l.J=l/2) 



- T 

3 D * N 



D*N 



+ 



„(I=l,J=3/2)' 



(5) 



where the T-matrices are evaluated at s = P - P , being Pq ~ qa + £jv(p) an< i 
P = q + p are the total energy and momentum of the meson-nucleon pair in the 
nuclear matter rest frame, and (qo,c[ ) and (Ej\f,p ) stand for the energy and momentum 
of the meson and nucleon, respectively, in this frame. The self-energy is determined 
self-consistently since it is obtained from the in-medium amplitude which contains the 
meson-baryon loop function, and this quantity itself is a function of the self-energy. 
Then, the meson spectral functions for D and D* mesons read 



Im77 D(£) .)(go,q) 



- m D(D*)- n D(D< 



)(<7o,q) 



(6) 



The D and D* spectral functions are displayed on the l.h.s. of Fig. [T] Those spectral 
functions show a rich spectrum of resonance-hole states. On one hand, the D meson 
quasiparticle peak mixes strongly with X' c (2823)iV~ 1 and X' C (2868)7V^ 1 states. On the 
other hand, the /lc(2595)iV -1 is clearly visible in the low-energy tail. With regard to 
the D* meson, the D* spectral function incorporates the J = 3/2 resonances, and the 
quasiparticle peak fully mixes with the E C (2902)N^ 1 and J 4 C (2941)/V" 1 states. For 
both mesons, the YciV" 1 modes tend to smear out and the spectral functions broaden 
with increasing phase space, as seen before in the SU(4) model [8]. Note, however, 
that not all the states described in Sec. [3] are seen in the D and D* spectral functions. 
This is due to the fact that some of those states do not strongly couple to DN and/or 
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Fig. 2 D°-nucleus bound states for different angular momentum L. 
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Fig. 3 D and D°- nucleus bound states for different angular momentum L. 



D* N states. Moreover, as we have just seen, resonances with higher masses than those 
described in Sec. [3] are also present in the spectral functions. Those resonant states 
were seen in the wider energy range explored in Ref . [16] . 



5 Charmed mesons in nuclei 

D and l)-meson bound states in 208 Pb were predicted in Ref. [26] relying upon an 
attractive D and D -meson potential in the nuclear medium. This potential was ob- 
tained within a quark-meson coupling (QMC) model [27] • The experimental observa- 
tion of those bound states, though, might be problematic since, even if there are bound 
states, their widths could be very large compared to the separation of the levels. This is 
indeed the case for the potential derived from a SU{4) t-vector meson exchange model 
for D-mesons [28] . 

We solve the Schrodinger equation in the local density approximation so as to 
analyze the formation of bound states with charmed mesons in nucleus. We use the 
energy dependent optical potential 

V(r,E)= n ^ = m+ ^ = °' (7) 

where E = q° — m is the D or D energy excluding its mass, and 77 the meson self- 
energy. The optical potential for different densities is displayed on the r.h.s of Fig. [1] 
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For D mesons we observe a strong energy dependence of the potential close to the D 
meson mass due to the mixing of the quasiparticle peak with the X' c (2823)iV™ 1 and 
X' c (2868)Af~ 1 states. As for the D meson, the presence of a bound state at 2805 MeV 
[17j . almost at DN threshold, makes the potential also strongly energy dependent. This 
is in sharp contrast to the SU(4) model (see analysis in Ref. [29]). 

Then, the question is whether D and/or D will be bound in nuclei. We start 
by discussing D mesons in nuclei. We observe that the D°-nucleus states are weakly 
bound (see Fig. [2]), in contrast to previous results using the QMC model. Moreover, 
those states have significant widths |29| . in particular, for 208 Pb |26| . Only D°-nucleus 
bound states are possible since the Coulomb interaction prevents the formation of 
observable bound states for D + mesons. 

Apropos of £>-mesic nuclei, not only D~ but also D° bind in nuclei as seen in 
Fig. [3] The spectrum contains states of atomic and of nuclear types for all nuclei for 
D~ while, as expected, only nuclear states are present for D° in nuclei. Compared to 
the pure Coulomb levels, the atomic states are less bound. The nuclear ones are more 
bound and may present a sizable width |30| . Moreover, nuclear states only exist for 
low angular momenta. 

The information on bound states is very valuable to gain some knowledge on the 
charmed meson-nucleus interaction, which is of interest for PANDA at FAIR. The 
experimental detection of D and D-meson bound states is, though, a difficult task. For 
example, reactions with antiprotons on nuclei for obtaining D°-nucleus states might 
have a very low production rate (see Ref. [29] for details). Reactions but with proton 
beams, although difficult, seem more likely to trap a D° in nuclei [29] . 



6 Conclusions and Future 

The properties of charmed mesons in infinite dense matter and nuclei have been studied 
within a unitary meson-baryon coupled-channel model which incorporates heavy-quark 
spin symmetry. Several resonances in the C = 1 and S — sector have been analyzed 
and compared with experimental data from several facilities as well as with other 
theoretical models. The in-medium modifications on those baryonic resonances and 
on charmed mesons have been obtained. The in-medium solution accounts for Pauli 
blocking effects and meson self-energies. We have analyzed the evolution with density 
of the open-charm meson spectral functions and finally studied the possible formation 
of D-mesic nuclei. On one hand, only weakly bound D°-nucleus states seem to be 
feasible. On the other hand, D~ and D°- nuclear bound states are possible, the latter 
ones with a stronger binding than for D° and with also a sizable width. The challenge 
is, though, the experimental detection. 
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